Electromagnetically Induced Transparency with Quantized Fields in Optocavity 
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We report electromagnetically induced transparency using quantized fields in optomechanical 
systems. The weak probe field is a narrow band squeezed field. We present a homodyne detection of 
EIT in the output quantum field. We find that the EIT dip exists even though the photon number 
in the squeezed vacuum is at the single photon level. The EIT with quantized fields can be seen 
even at temperatures of the order of 100 mK paving the way for using optomechanical systems as 
memory elements. 
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I. INTRODUCTION 

The interaction of a nano mechanical system via the 
radiation pressure [3, 0| is like a three wave interaction 
in nonlinear optics l3l|. This interaction can lead to pro- 
cesses like upconversion for example a photon of fre- 
quency LUc can be converted into a photon of frequency 
LOp = uJc + ojm where LUm is the frequency of the mechan- 
ical oscillator. Such upconversion processes have been 
useful in cooling the nano mechanical systems 0-01 ■ In 
a previous paper ^ we showed how such upconversion 
processes can lead to electromagnetically induced trans- 
parency in optomechanical systems. The EIT in such 
systems turned out to share many of the features of EIT 
in atomic vapors. The EIT in o pto mechanical systems 
has been seen experimentally [Mt11|- Traditionally al- 
most all EIT experiments in atomic systems and other 
systems have been done with coherent pump and probe 
fields [T2| - [l^ . Akamatsu et al [11] did the very first ex- 
periment on EIT using the squeezed light in atomic va- 
pors. They essentially reported that squeezing of the 
probe is not graded much by the quantum noise of the 
medium under EIT conditions. Following this a num- 
ber of other experiments [Tg', on EIT using quantized 
fields were reported. The EIT with quantized fields is 
very significant in storage of fields at single photon level 

MM. 



In this paper we examine EIT in optomechanical sys- 
tems using the quantized fields. In optomechanical sys- 
tems the noise is added both by the resonator's noise as 
well as the noise of the mechanical system. We would 
find conditions when perfect EIT of the quantized field 
results. We would study how the temperature of the me- 
chanical system can degrade the EIT. We would present 
detailed results for designs of nano mechanical systems as 
used in refs [1, [l^ . We find that certain designs of nano 
mechanical systems are good even at temperatures of the 
order of lOOmK. Thus such systems would be quite useful 
as optical memories at single photon level. The results 
that we present can be extended to reactive case [23l - l25| . 

The organization of the paper is as follows. In Sec. II, 
we describe the model, derive the equations of the motion 



for the system, and obtain the steady-state mean values. 
In Sec. Ill, we show how to detect EIT with quantized 
fields, we present a homodyne detection and obtain the 
relevant spectrum. In Sec. IV, we discuss the impact 
of the coupling field on the homodyne spectrum of the 
output field, and show the existence of the EIT in the 
homodyne spectrum of the quantized field at the output. 
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FIG. 1: Sketch of the studied system. A coherent coupling 
field at frequency uJc and a squeezed vacuum at frequency cup 
enter the cavity through the partially transmitting mirror. 



II. MODEL 



The model we are going to consider has been discussed 



in detail before [26|, |27| and is sketched in Fig. 1. The 

cavity consists of a fixed mirror and a movable mirror 
separated by a distance L. The fixed mirror is partial 
transmitting, while the movable mirror is 100% reflect- 
ing. The cavity is driven by a strong coupling field at 
frequency uJc- A quantized weak probe field in a squeezed 
vacuum state at frequency ujp is injected into the cavity 
through the fixed mirror. The movable mirror interacts 
with the cavity field through the radiation pressure. The 
movable mirror is modeled as a harmonic oscillator with 
mass TO, frequency tUm, and decay rate 7™- Moreover, 
the movable mirror and its environment are in thermal 
equilibrium at a low temperature T. 
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In such a system, the couphng between the movable 
mirror and the cavity field is dispersive so that the fre- 
quency LOoiq) of the cavity field depends on the displace- 
ment q of the movable mirror: ujQ{q) — mTc/{L+q), where 
c is the light speed in vacuum, and n is the mode number 
in the cavity. For q ^ L, we can expand u^oiq) to the 
first order of q, thus we have ujoiq) ~ t^o(O) + ^^^§^9 ~ 
Wo — ^q, where we write a;o(0) as luq. 

Let c (c^) be the annihilation (creation) operators for 
the cavity field, Q and P be the dimensionless operators 
for the position and momentum of the movable mirror 

with Q = ^Ini^q and P \J^^-P- Note that the 

commutation relation for Q and P is [Q, P] = 2i. In a 
frame rotating at the frequency ujc of the coupling field, 
the Hamiltonian for the system is 



H = h{LUQ — LUc)c^c — hgc) cQ 
-\-ihe{c^ — c). 



hujri 



■(g' + p") 



(1) 



In the above equation, the parameter g — 
{ujc/ L) h/ {2muj,n) is the coupling strength between the 
cavity field and the movable mirror, where we assume 
Wo — The parameter e is the real amplitude of the 

coupling field, depending on its power p by e = \Jj^7 
where k is the photon loss rate due to the transmission 
of the fixed mirror. 

The time evolution of the total system is obtained from 
Hamiltonian ([T]) by deriving the Heisenberg equations of 
motion and adding the damping and noise terms. The 
basic equations are given by 



Q ^ WmP, 

P = 2gnc - UJmQ - IrnP + 

c = i{ujc — tjJo + gQ)c + e — KC+ yplkcim 
ct = — i(wc — Wo + gQ)c^ + e — Kct + \/2k4„ 



(2) 



Here we have introduced the thermal Langevin force ^ 
with vanishing mean value, resulting from the coupling 
of the movable mirror to the environment. The Langevin 
force ^ has the correlation function in the frequency do- 
main: 



(e(w)e(l^)) = 47r7„ 



w 



1 -f coth 



/ hijj 



\2KbT 



(3) 



where Kb is the Boltzmann constant. Throughout this 
paper the following Fourier relations are used 



1 ;'+oc 

fit) - ^ / f{^)e-^^'du, 

fit) - ^ / fi 



w)e-*'"*dw. 



(4) 



where f^—ui) = [f{—uj)]^. The represents the input 
quantum field which is centered around the frequency 
LUp ~ ujc + LOm with a finite bandwidth F. The quantized 
field has the following nonvanishing correlation functions. 



F2 -f (w - W™)2 

r2 + (w - W™)2 



(5(w -I- — 2Wm), 



(5(w + f}), 



(5) 



where N is the photon number in the squeezed vac- 
uum, and M = y/N{N + 1). The antinormally ordered 
term has a broad band contribution coming from vac- 
uum noise. Note that by setting M — 0, we would 
obtain standard phase independent quantum field with 
mean number of photons r^+(^-Lj )^ around the fre- 
quency w = W„j . 

The mean values at steady state can be obtained from 
Eq. ^ by setting all time derivatives to zero. These are 
found to be 



Ps = 0, Qs 



2g\cs 



UJ„ 



where 



A = Wo - Wc - gQs 
is the effective cavity detuning. 



(6) 



(7) 



III. 



THE OUTPUT FIELD AND ITS 
MEASUREMENT 




FIG. 2: Sketch of the measurement of the output field. The 
output field Coutit) is mixed with a strong local field CLoit) 
centered around the probe frequency at a beam splitter, 
where Coutit) is defined as the sum of the output field Cout(t) 
from the cavity and the input quantized field Cin{t). BS: 50/50 
beam splitter, PD: photodetector, SA: spectrum analyzer. 

The output field is a quantum field, it contains many 
Fourier components. Since the quantized input field is 
centered around Wp = Wc-I-Wm, the interesting component 
of the output field is near the probe frequency w„. So we 
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mix the output field Cout{t) with a strong local field CLo{t) 
centered around the probe frequency iOp at a 50/50 beam 
splitter, as shown in Fig. 2. In a frame rotating at the 
frequency oJc, CLo{t) = CLoe~^^°*, where Sq = ojp — ujc- 
The difference between the output signals from the two 
photodetectors is sent to the spectrum analyzer, and the 
output signal from the spectrum analyzer depends on the 
phase of clo- If clo is real, the homodyne spectrum ^(a;) 
of the output field measured by the spectrum analyzer is 
given by 

{[c*Loit)£outit) + C.C.][cl,it')dout{t') + C.C.]) 

- ^ J dcje-*"(*-*')X(cj). (8) 



Thus in our investigations of EIT with quantized fields, 
X{uj) is the quantity of interest. 

In order to study the EIT effect in the homodyne spec- 
trum X{ijj) of the output field, we will calculate the fiuc- 
tuations of the output field. The steady state part would 
not contribute as it is at the frequency of the coupling 
field. We assume the photon number in the cavity is 
large enough so that each operator can be written as a 
linear sum of the steady-state mean value and a small 
fluctuation, which yields 

Q = Qs+SQ, P = Ps+6P, c = Cs+Sc, (9) 

where 5Q, SP, and Sc are the small fluctuations around 
the steady state. By substituting Eq. ([9]) into Eq. 

one can arrive at the linearized equations for the 
fluctuation operators. Further we transform the lin- 
earized equations into the frequency domain by Eq. ^ 
and solve it, we can obtain the fluctuations dc{uj) of 
the cavity field. Then using the input-output relation 
Cout{i^) = V2kc{(jj) — Cin{uj), wc Can find the fluctuations 
(5cout(w) of the output held. For the purpose of Fig. ^ 
we deflne the output field as Couti>jj) — Cout(w) -I- Cm(w), 
then we find the result 

= Viu)aLo) + E{uj)c,r.{uj) + F{u)cl{-uj), (10) 



2k 

x[k- i{uj + A)]} 

4k 2 2 

——UJrng C^l, 

d(uj) 



in which 

^(^) = 37-T'^™ff^c^i, (11) 
where 



d{uj) = -4a;,„A52|c3p + (w^-cj2-i7„ja;)[(K-ia;)^-f A^]. 

(12) 

The first term on the right-hand side of Eq. ((TU)) refers 
to the contribution of the thermal noise of the movable 
mirror, and the other two terms represent the contribu- 
tion of the squeezed vacuum. To illustrate the mean- 
ing of the last two terms, let the squeezed vacuum be 
a single mode i.e. Cmit) — (7e~'('^p~'^'=)*, then Ci„(aj) — 
2ttCS{ll! — ujp + ujc) and cj„(— w) = 2ttC''6{uj + ojp — ojc). 
Thus the fluctuations of the output field 8cout{t) — 
^j+^V{Lo)i{u:)e-^^'duj + CE{ujp - u;,)e-'(-''-"=)* + 

CtF(wc - cjp)e~*("'=^"p)*. Therefore, E{ujp - ujc) is the 
component at the probe frequency ujp [which in rotating 
frame is Wp — wj, and F{uJc — ^p) is the component at 
the new frequency 2ujc ~ cOp [which in rotating frame is 

UJ, 

movable mirror and the cavity field 



ujp], due to the nonlinear interaction between the 



By the aid of the correlation functions of the noise 
operators Ci„(a;) and C(w), and neglecting fast oscillat- 
ing terms at frequency ±2cjm, we obtain the homodyne 
spectrum X{lu) of the output field as measured by the 
set up of Fig. [2] 



X(w) 



E{uJ+UJ„^)E{-UJ+UJ.,n)z-:^, K + \E{uJ + OJm)\^ 



+ \E{-OJ + OJm)\ 



+E [-UJ+UJrn)E [UJ + UJ„i)—^ ' ' 



+ \E{lU + UJ,n)r + \F{-UJ + UJ„,W 

-fin^+^„)|227™-^^'^" 

+ \V{~UJ + iJ,n)f2-f„ 



UJ — LOr, 



|l + coth 


'h{uj + UJrnY 




2kBT 



m 



1 + coth 



h{u - uj„i) 

2kBT 



(13) 



where the first four terms in Eq. (1131) originate from the squeezed vacuum, the next two terms not involving N 
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and M are the contributions of the spontaneous emission 
of the input vacuum noise; the last two terms in Eq. (jl3p 
result from the thermal noise of the movable mirror. 

IV. EIT IN THE HOMODYNE SPECTRUM OF 
THE OUTPUT QUANTIZED FIELD 

After having derived the homodyne spectrum of the 
output field, we next examine it numerically to explore 



We use the parameters from the experimental paper [9| 
focusing on the EIT in the optomechanical system: the 
wavelength of the coupling field A = ^-ncjujc — 775 nm, 
the coupling constant g = 27r x 12 GHz/nmy^fiTI^^ij^, 
the mass of the movable mirror m = 20 ng, the frequency 
of the movable mirror — 27rx 51.8 MHz, the cavity de- 
cay rate k = 27rx 15 MHz, njujm — 0.289, the mechanical 
damping rate 7^ = 27r x 41 kHz, the mechanical quality 
factor Q' = uJmhm = 1263. In addition, we choose the 
linewidth of the squeezed vacuum F = 2k, and consider 
the resonant case A = Wm- 

For = 10 and M = y/N{N +1) « 10, p = 20 mW, 
T = 20 mK, the first term in Eq. (fT4| which is the con- 
tribution of the squeezed vacuum is about 6.5 x 10~*, the 
sum of the second and third terms in Eq. (fT4|) which are 
the contribution of the input vacuum noise is about 0.16, 
the last term arising from the thermal noise of the mov- 
able mirror is about 0.14. The contribution of the input 
quantum field in principle can be obtained by doing the 
experiment with and without the quantized field and by 
subtracting the data i.e. by studying X{Q) — X{Q)\n=o. 
The squeezed field part in a sense exhibits perfect EIT. 
If M = 0, i.e. the input quantized field is phase insensi- 
tive, then such a field leads to a term 2N\E{u)m)\'^ which 
is equal to 1.6 for the above mentioned parameters and 
hence there is no perfect EIT. The squeezed field changes 
2N\E{uj^)\^ to N[E{u]^) + E* {oj„,)? and for the above 
parameters the number changes from 1.6 to 6.5 x 10~^. 

For = 5, M = y/N{N + 1) and 0, T = 20 mK, we 
plot the homodyne spectrum X{u)) of the output field as 
a function of the normalized frequency oj /iOm in the ab- 
sence (dotted curve) and presence (solid, dotdashed, and 
dashed curves) of the coupling field in Fig. |3l First let us 
look at the case that the input quantum field is phase de- 
pendent {M ~ y/N{N + 1)). In the absence of the cou- 
pling field, one can note that the homodyne spectrum of 
the output field has a Lorentzian line shape. However, in 
the presence of the coupling field at different power levels, 
the solid curve (p = 10 mW and M = y^N{N + 1)) and 



the EIT phenomenon in the homodyne spectrum of the 
output field. Since the original equations ([2]) are nonlin- 
ear, these can have instabilities. Thus in the following, 
we work in the stable regime of the system. We first 
examine the frequency at which we expect transparency. 
This is w = 0. For A^ « M, 




FIG. 3: Plots of the homodyne spectrum X{ljj) as a function 
of Lj/ojm for A'^ = 5 in the absence (dotted curve) and the 
presence (solid, dotdashed, and dashed curves) of the coupling 
field for the temperature of the environment T — 20 mK. The 
solid curve is for p = 10 mW and M = \/N{N + 1), the 
dotdashed curve is for p = 20 mW and M = -yAXiVTTy, 
the dashed curve is for p = 20 mW and M — 0. 



the dotdashed curve (p = 20 mW and M = y/N{N + 1)) 
exhibit the EIT dip, which is the result of the destruc- 
tive interference between the squeezed vacuum and the 
scattering quantum field at the probe frequency ujp gen- 
erated by the interaction of the coupling field with the 
movable mirror. For p = 20 mW and M = y/ N{N + 1), 
the minimum value of X{uj) is about 0.22. Moreover, the 
linewidth of the dip for p = 20 mW is larger than that 
for p = 10 mW due to power broadening. Generally the 
EIT dip has a contribution to its width which is propor- 
tional to the power of the coupling field. We indeed find 
that the width for p = 20 mW is 0.26wm, which is about 
twice the width for p = 10 mW. If the input quantum 
field is phase independent (M = 0) [the dashed curve], 
we can see that the maximum value of X{uj) for p = 20 
mW and M = is about half that for p = 20 mW and 



X{0) - N[E{L0,n)+E^UJ,n)? + \E{uj,n)\^ + \F{uj^)\^+A\VM\^J,nCOth 
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M = y/N{N + 1). 

Next we increase the temperature to 100 mK. The fig- 
ure [4] displays the homodyne spectrum X{uj) of the out- 
put field against the normalized frequency Lo/ujm in the 
absence (dotted curves) and presence (solid curves) of 
the coupling field for = 1,5 and M = ^N(N + 1). 
In the presence of the coupling field (p = 10 mW), it is 
seen that the EIT dip still appears in the homodyne spec- 
trum of the output field for = 1 and 5. Note that the 
two dips almost have the same minimum values (about 
1.43) and the same linewidth (about 0.15w,„). Hence the 
temperature of the environment is not detrimental to the 
EIT behavior. 



X(cj) 




FIG. 4: Plots of the homodyne spectrum X{lo) as a function 
of ui/uJm for different values of the parameter A'^ and M = 
^JN{N + 1) in the absence (dotted curves) and the presence 
(solid curves) of the coupling field with power p = 10 mW 
and the temperature of the environment T = 100 mK. The 
upper two curves are for = 5, the lower two curves are for 
N = 1. 

The effects discussed above occur under wide range of 
parameters. We demonstrate this by using the experi- 
mental parameters [11]: A — 27rc/wc — 1064 nm, L = 25 
mm, 5 « 27r X 11.28 MHz/nm-^/ h/ {2mu}„i) , m = 145 ng, 
u},n = 27r X 947 kHz, k = 27r x 215 kHz, = 0.227, 

7m = 27r X 141 Hz, Q' = LOm/'Jm = 6700. The values for 
parameters T, p, A^, M , F, and A are the same as those in 
Fig. 21 Shown in Fig. [S]is the homodyne spectrum X{lu) 
of the output field as the normalized frequency uj /ujm is 
varied for T ^ 100 mK and p 0, 10 mW. Note that the 



EIT exists for A" = 1 and 5 in the presence of the cou- 
pling field. The linewidth of the dip for A^ = 5 is about 
0.2u}m, and as expected gets broadened due to power. We 
have further studied the effect of temperature and we find 
that there is rather weak dependence of the EIT curves 
on temperature. Therefore, current optomechanical de- 
signs can be used to realize quantum optical memory at 
single photon level. 

X(6J) 




-0.4 -0.2 0.2 0.4 

FIG. 5: Plots of the homodyne spectrum X{ijj) as a function 
of LojiOra for different values of the parameter X and M = 
\J X(X -I- 1) in the absence (dotted curves) and the presence 
(solid curves) of the coupling field with power p = 10 mW 
and the temperature of the environment T = 100 mK. The 
upper two curves are for A' = 5, the lower two curves are for 
X 



V. CONCLUSIONS 

In conclusion, we have demonstrated EIT using quan- 
tum fields in optomechanical systems under a wide range 
of conditions. For squeezed quantum fields we obtained 
perfect EIT. The EIT gets degraded in phase insensitive 
quantum fields. We have shown that even temperature 
is not critical for observations of EIT. The results get be 
generalized to optomechanical systems working on the re- 
active coupling [23l - [25j . Our work suggest that optome- 
chanical systems could be used as elements for quantum 
memory. 

We gratefully acknowledge support from the NSF 
Grant No. PHYS 0653494. 



[1] P. Meystre, E. M. Wright, J. D. McCuUen, and E. Vignes, 

J. Opt. Soc. Am. B 2, 1830 (1985); M. Bhattacharya, H. 

Uys, and P. Meystre, Phys. Rev. A 77, 033819 (2008). 
[2] S. Mancini, V. I. Manko, and P. Tombesi, Phys. Rev. A 

55, 3042 (1997). 
[3] R. W. Boyd, Xonlinear Optics, 2nd ed. (Academic, San 

Diego, 2003). 

[4] S. Gigan, H. R. Bohm, M. Paternostro, F. Blaser, G. 
Langer, J. B.Hertzberg, K. C. Schwab, D. Bauerle, M. 



Aspelmeyer, and A. Zeilinger, Nature (London) 444, 67 
(2006). 

[5] O. Arcizet, P. -F. Cohadon, T. Briant, M. Pinard, and 
A. Heidmann, Nature (London) 444, 71 (2006). 

[6] A. Schliesser, R. Riviere, G. Anetsberger, O. Arcizet, and 
T.J. Kippenberg, Nature Physics 4, 415 (2008). 

[7] Y. Park and H. Wang, Nature Physics 5, 489 (2009). 

[8] G. S. Agarwal and S. Huang, Phys. Rev. A 81, 041803(R) 
(2010). 



6 



[9] S. Weis, R. Riviere. S. Deleglise, E. Gavartin, O. Arcizet, 
A. Schliesser, and T. J. Kippenberg, Science 300, 1520 
(2010). 

[10] Q. Lin, J. Rosenberg, D. Chang, R. Camacho, M. Eichen- 
ficld, K. J. Vahala, and O. Painter, Nature Photonics 4, 
236 (2010). 

[11] A. H. Safavi-Naeini, T. P. Mayer Alegre, J. Chan, M. 

Eichenfield, M. Winger, Q. Lin, J. T. Hill, D. Chang, O. 

Painter, arXiv:1012.1934 
[12] S. E. Harris, J. E. Field, and A. Imamoglu, Phys. Rev. 

Lett. 64, 1107 (1990). 
[13] K. -J. Boiler, A. Imamoglu, and S. E. Harris, Phys. Rev. 

Lett. 66, 2593 (1991). 
[14] O. Kocharovskaya, Y. Rostovtsev, and M. O. Scully, 

Phys. Rev. Lett. 86, 628 (2001). 
[15] D. Akamatsu, K. Akiba, and M. Kozuma, Phys. Rev. 

Lett. 92, 203602 (2004). 
[16] D. Akamatsu, Y. Yokoi, M. Arikawa, S. Nagatsuka, T. 

Tanimura, A. Furusawa, and M. Kozuma, Phys. Rev. 

Lett. 99, 153602 (2007). 
[17] M. Arikawa, K. Honda, D. Akamatsu, Y. Yokoil, K. Ak- 
iba, S. Nagatsuka, A. Furusawa, and M. Kozuma, Optics 



Express, 15, 11849 (2007). 
[18] M. D. Eisaman, A. Andre, F. Massou, M. Fleischhauer, 

A. S. Zibrov, and M. D. Lukin, Nature 438, 837 (2005). 
[19] J. Appel, E. Figueroa, D. Korystov, M. Lobino, and A. 

L Lvovsky, Phys. Rev. Lett. 100, 093602 (2008). 
[20] M. Lobino, C. Kupchak, E. Figueroa, and A. L Lvovsky, 

Phys. Rev. Lett. 102, 203601 (2009). 
[21] K. Honda, D. Akamatsu, M. Arikawa, Y. Yokoi, K. Ak- 
iba, S. Nagatsuka, T. Tanimura, A. Furusawa, and M. 

Kozuma, Phys. Rev. Lett. 100, 093601 (2008). 
[22] S. Groblacher, K. Hammerer, M. Vanner, and M. As- 

pelmeyer. Nature (London) 460, 724 (2009). 
[23] M. Li, W. H. P. Pernice, and H. X. Tang, Phys. Rev. 

Lett. 103, 223901 (2009). 
[24] S. Huang and G. S. Agarwal, Phys. Rev. A 81, 053810 

(2010). 

[25] S. Huang and G. S. Agarwal, Phys. Rev. A 82, 033811 
(2010). 

[26] K. Jahne, C. Genes, K. Hammerer, M. Wallquist, E. S. 

Polzik, and P. ZoUer, Phys. Rev. A 79, 063819 (2009). 
[27] S. Huang and G. S. Agarwal. arXiv:0905.4234 , 



